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Abstract.  In  1994  June-July  we  obtained  nearly  concur¬ 
rent  measurements  of  Cyg  X-3  in  hard  X-rays  and  the 
infrared,  using  the  OSSE  instrument  on  GRO  and  the 
UKIRT.  OSSE  made  a  high  significance  measurement  of 
the  hard  X-ray  (>  50  keV)  orbital  light  curve.  Five  days 
after  the  end  of  the  OSSE  observations  we  obtained  a 
good  IR  K-band  light  curve  from  data  covering  one  orbit 
of  Cyg  X-3.  Measurements  of  the  light  curve  phase  and 
shape  probe  the  geometry  of  the  system  and  the  physics 
of  radiation  production  and  scattering.  We  determined 
for  the  first  time  that  the  light  curve  phases  are  consis¬ 
tent  at  these  two  different  wavelengths.  The  light  curve 
shapes  were  both  marginally  inconsistent  (at  about  the 
0.2%  level)  with  the  canonical  soft  X-ray  shape  measured 
by  EXOSAT  in  1983-1985.  Underlying  variations  in  DC 
IR  flux  and  flaring  during  the  observation  make  it  diffi¬ 
cult  to  draw  unambiguous  conclusions;  more  observations, 
with  better  IR  coverage,  are  needed. 

Key  words:  X-ray  binaries  -  Cygnus  X-3  -  X-rays  - 
Infrared 


1.  Introduction 

One  of  the  mysteries  associated  with  the  bright  X-ray 
source  Cyg  X-3  is  the  change  in  the  ~  4.8  hr  period.  The 
period,  commonly  thought  to  correspond  to  the  orbital  cy¬ 
cle  of  a  compact  object /companion  star  system,  is  increas¬ 
ing,  with  P/P  ~  10-6  yr-1  (e.g.,  Kitamoto  et  al.  1992). 
This  is  inconsistent  with  the  Roche  lobe  transfer  of  mate¬ 
rial  from  a  main  sequence  star  to  the  compact  object  (Mol- 
nar  1988).  It  can  be  explained,  under  certain  conditions, 
by  a  wind  from  the  companion  star  (van  Kerkwijk  et  al. 
1992).  The  period  is  typical  of  low-mass  X-ray  binaries; 
however,  there  is  evidence  that  the  companion  is  actually 
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a  high-mass  Wolf-Rayet  star  (van  Kerkwijk  et  al.  1992). 
There  is  also  a  significant  (4<r)  P  term  in  the  ephemeris, 
first  measured  by  Van  der  Klis  and  Bonnet-Bidaud  (1989). 

The  orbital  light  curve,  observed  in  IR  and  X-rays,  is 
asymmetric.  This  requires  some  asymmetry  in  the  system, 
either  in  the  binary  orbit  itself  (i.e. ,  a  significant  eccen¬ 
tricity)  or  in  the  matter  in  the  system.  If  the  asymmetry 
is  caused  by  an  elliptical  orbit,  the  light  curve  shape  is 
expected  to  change  significantly  over  time  due  to  apsidal 
motion  (periastron  advance)  (Ghosh  et  al.  1981).  If  the  or¬ 
bit  is  elliptical,  apsidal  motion  may  also  be  responsible  for 
some  or  all  of  the  apparent  period  increase  (Eisner  et  al. 
1980).  Therefore,  measurements  of  the  light  curve  shape 
are  directly  related  to  the  orbital  dynamics.  Short-term 
variations  in  the  light  curve  shape,  however,  complicate 
the  interpretation  of  individual  observations  (e.g.,  van  der 
Klis  &:  Bonnet-Bidaud  1981). 

The  asymmetric,  non-eclipsing  light  curve,  along  with 
the  observed  He  lines  in  the  IR  spectrum  (van  Kerkwijk 
et  al.  1992)  imply  that  the  primary  source  is  surrounded 
by  dense  material.  This  material  may  be  in  the  form  of  a 
wind  from  the  secondary  (Pringle  1974;  Willingale,  King, 
&:  Pounds  1985),  a  “cocoon”  (Milgrom  &:  Pines  1978),  or 
an  accretion  disk  corona  (White  &:  Holt  1982).  See  Bonnet- 
Bidaud  &:  Chardin  (1988)  for  a  valuable  review  of  many 
of  the  Cyg  X-3  observations  and  theories. 

It  has  been  theorized  that  Cyg  X-3  shows  so  many 
unique  properties  because  it  is  in  a  very  short-lived,  tran¬ 
sitional  phase  of  binary  evolution.  If  this  is  the  case  a 
study  of  its  dynamics  (via  the  orbital  ephemeris),  and  lo¬ 
cal  environment  (via  the  light  curve  and  spectrum)  may 
have  broad  implications  for  X-ray  binaries  in  general. 

2.  Observation  Summary 

2.1.  OSSE 

The  OSSE  instrument  on  GRO  is  described  in  detail  by 
Johnson  et  al.  (1993).  It  consists  of  four  independent  cylin¬ 
drical  Nal/Csl  phoswhich  scintillators,  actively  shielded 
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on  the  sides  by  Nal.  The  normal  spectral  energy  range 
is  50  keV  to  10  MeV.  The  front  of  each  detector  is  cov¬ 
ered  by  a  tungsten  collimator,  which  defines  a  rectangular 
(11.4°  x  3.8°  FWHM)  field  of  view.  Each  detector  may 
be  positioned  independently  about  an  axis  parallel  to  the 
long  direction  of  the  collimator.  Spectral  data  are  accumu¬ 
lated  in  two-minute  intervals,  alternating  between  source 
and  background  pointings  of  the  detectors.  In  the  Cyg  X- 
3  observations  the  orientation  of  the  spacecraft  and  the 
positions  of  the  detector  pointings  were  chosen  so  that 
there  was  no  significant  flux  contribution  from  Cyg  X-l 
in  either  the  source  or  the  background. 

OSSE  made  seven  significant  observations  of  Cyg  X-3 
from  1991  to  1994:  1991  May  30  -  Jun  15  (viewing  period 
2),  1991  Aug  8-15  (VP  7),  1991  Nov  28  -  Dec  11  (VP  15), 
1992  Dec  9-15  (VP  203),  1994  May  24  -  31  (VP  328), 
1994  June  7-17  (VP  330-331),  and  1994  July  7-12  (VP 
333).  OSSE  has  detected  significant  flux  from  Cyg  X-3 
up  to  ~  200  keV.  Some  results  of  the  earlier  observations 
have  been  described  by  Matz  et  al.  (1994).  We  focus  here 
on  the  measurements  of  Cyg  X-3  taken  during  1994  June 
7-17  and  July  7-12  (TJD  9510-9520  and  9540-9545). 

2.2.  UKIRT 

Cyg  X-3  was  observed  with  the  UK  Infrared  Telescope 
(UKIRT)  on  the  night  of  1994  July  17  UT  (TJD  9550) 
over  nearly  one  entire  4.8  hr  orbital  cycle,  using  the  IR- 
CAM3  instrument.  IRCAM3  is  a  256  x  256  imaging  ar¬ 
ray,  operating  in  the  1-5  micron  range,  with  a  resolution 
of  0.286  arcsec/pixel.  The  source  was  monitored  in  the 
J  (1.2  micron),  H  (1.6  micron),  K  (2.2  micron)  and  nbL 
(3.6  micron)  bands.  Weather  conditions  were  not  good  but 
seeing  was  reasonable  and  nearby  field  stars  allowed  accu¬ 
rate  relative  photometry  to  be  performed.  Faint  standard 
G93-48  was  used  for  absolute  flux  calibration  at  J,  H,  and 
K.  A  comparison  of  variability  in  all  four  bands  will  be 
discussed  elsewhere.  Mosaicing  of  images  revealed  at  least 
ten  previously  undetected  objects  within  a  few  arcsec  of 
Cyg  X-3  —  see  Fender  &:  Bell  Burnell  (1995)  for  a  new 
K-band  finder  chart. 

3.  OSSE  Ephemeris  Derivation 

The  data  in  this  paper  were  analysed  using  an  up-to-date 
X-ray  ephemeris  derived  from  previous  observations  (see 
Kitamoto  et  al.  (1992)  and  references  therein),  the  last  of 
which  was  in  1987,  plus  the  seven  OSSE  observations  from 
1991  to  1994. 

For  each  of  the  seven  OSSE  observations  of  Cyg  X-3 
the  two-minute  background-subtracted  count  rate  spectra 
were  integrated  over  the  energy  range  44-130  keV.  The 
time  of  each  two-minute  observation  was  corrected  to  the 
solar  system  barycenter.  The  orbital  phase  was  then  cal¬ 
culated;  as  a  test  of  ephemeris  dependence,  three  different 
methods  were  used:  a  constant  period  (0.19969079d)  and 


Time  (hrs  UT) 

Fig.  1.  UKIRT  corrected  K-band  data  from  July  17,  1994. 
Individual  points  represent  10  s  accumulations.  Error  bars  are 
roughly  consistent  with  the  symbol  size. 

the  quadratic  and  cubic  ephemerides  of  van  der  Klis  &: 
Bonnet- Bidaud  (1989).  Orbital  light  curves  were  produced 
by  accumulating  the  data  into  the  appropriate  phase  bins. 
Due  to  the  relatively  weak  flux  from  Cyg  X-3  at  hard  X- 
ray  energies,  the  OSSE  data  must  be  epoch-folded  over 
a  number  of  days  to  produce  a  significant  signal.  This 
has  the  advantage  of  averaging  over  the  known  cycle-to- 
cycle  variations  in  the  light  curve  shape.  A  significant 
orbital  modulation  of  the  flux  was  observed  in  each  of 
the  seven  viewing  periods  except  VP  203,  where  the  fit¬ 
ted  light  curve  amplitude  was  only  marginally  positive 
(~  2.8(7  above  zero). 

The  EXOSAT  template  light  curve  (van  der  Klis  &: 
Bonnet- Bidaud  1989)  was  fit  to  the  OSSE  light  curve 
for  each  observation,  allowing  the  phase,  amplitude,  and 
background  (DC)  level  to  vary.  In  all  cases  except  for  VP 
333  the  X-ray  light  curve  template  provided  a  statistically 
acceptable  fit  to  the  OSSE  data;  in  VP  333  the  %2  was  high 
(58  for  29  dof)  but  the  template  seemed  to  describe  the 
data  well  enough  to  determine  the  phase.  The  template  is 
better  than  the  alternative  model,  a  simple  sine  function, 
which  has  commonly  been  used  to  fit  the  X-ray  data  to 
determine  pulse  arrival  times  even  though  it  is  not  a  good 
description  of  the  average  X-ray  light  curve  shape. 

For  each  of  the  light  curves,  the  shift  of  the  fitted 
template  from  phase  0  as  predicted  by  the  ephemerides 
was  used  to  calculate  a  pulse  arrival  time  for  a  cycle  near 
the  center  of  the  observation.  The  pulse  arrival  time  is 
equal  to  the  predicted  arrival  time  corrected  by  the  fit¬ 
ted  difference  in  phase  multiplied  by  the  period.  There 
are  no  significant  differences  between  the  pulse  arrival 
times  derived  using  the  quadratic,  cubic,  and  constant  pe¬ 
riod  ephemerides.  The  OSSE  arrival  times  are  consistent 
(within  errors)  with  extrapolations  of  the  ephemerides  of 
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Kitamoto  et  al.  (1992)  and  van  der  Klis  &:  Bonnet-Bidaud 
(1989). 

We  then  derived  a  new  ephemeris  by  combining  the 
OSSE  arrival  times  with  the  earlier  soft  X-ray  results. 
Consistent  with  previous  fits  (e.g.,  van  der  Klis  &:  Bonnet- 
Bidaud  1989;  Kitamoto  et  al.  1992),  we  can  exclude  a  sim¬ 
ple  quadratic  ephemeris  (at  a  confidence  level  of  10“ 11  or 
better),  while  a  cubic  description  (with  a  non-zero  P)  does 
fit  the  data. 

The  preliminary  results  are  shown  in  Table  1,  with 
parameters  defined  as  in  van  der  Klis  &:  Bonnet-Bidaud 
(1989).  Tn  is  the  arrival  time  of  the  minimum  of  the  nth 
orbital  cycle.  The  numbers  in  the  final  OSSE  ephemeris 
may  change  slightly  with  refinements  in  the  phase  deter¬ 
mination. 

Table  1.  PRELIMINARY  Best-Fit  Cubic  Ephemeris 
(T„  =  T0  +  Pon  ±  c0n 2  +  dn 3) 


0.0  0.5  1.0  1.5  2.0 

Phase 

Fig.  2.  OSSE  and  UKIRT  light  curves  with  best-fit  X-ray 
template  overplotted  (dashed  line).  Two  cycles  are  shown. 


Parameter 

Fitted  value 

To 

2440949.8999  ±  0.0014  HJD 

To 

0.19968240  ±  (2.8  x  10“7)  d 

Co 

(1.81  ±0.16)  x  10“10  d 

d 

(-2.04  ±0.25)  x  10“15  d 

X2  (dof) 

114.68  (84) 

Conf.  level 

0.015 

4.  Orbital  light  curve  analysis 

To  produce  the  orbital  IR  and  X-ray  light  curves  for  the 
1994  June-July  observations,  the  two  flaring  periods  (Fig. 
1)  were  first  removed  from  the  IR  data.  The  data  were  then 
corrected  to  the  solar  system  barycenter  (SSB)  and  epoch- 
folded  using  orbital  phases  calculated  from  the  OSSE  cu¬ 
bic  ephemeris  shown  in  Table  1.  For  comparison,  the  IR 
and  OSSE  data  were  binned  at  same  time  resolution.  The 
plotted  errors  in  the  IR  data  are  not  the  statistical  errors 
of  the  individual  measurements,  but  the  observed  variance 
of  the  data  included  in  each  phase  bin.  The  size  of  these 
variances  reflect  real  source  fluctuations,  as  determined 
from  relative  photometry  referenced  to  nearby  stars  Z  and 
D. 

The  resulting  light  curves  were  fitted  with  the  EX- 
OSAT  soft  X-ray  template  (van  der  Klis  &:  Bonnet-Bidaud 
1989)  in  the  manner  described  above.  The  phase  and  am¬ 
plitude  of  the  template,  and  the  DC  background  intensity 
were  varied  to  minimize  y2  for  the  X-ray  and  IR  data 
independently.  Figure  2  shows  the  resulting  light  curves 
with  best-fit  X-ray  template  (dashed  line)  overplotted. 
The  gaps  in  the  IR  light  curve  correspond  to  gaps  in  the 
data  collection.  From  the  fits  we  can  determine,  for  each 
of  the  two  data  sets,  the  phase  of  the  minimum  and  the 


consistency  of  the  measured  shape  with  the  soft  X-ray 
template. 

5.  Results 

5.1.  Light  curve  phase 

The  fitted  OSSE  phase  is  0.016±0.022  and  the  IR  phase  is 
0.069  ±0.024.  The  hard  X-ray  and  IR  light  curves  are  con¬ 
sistent  with  the  same  phase  at  minimum  within  1.6(7.  Ma¬ 
son,  Cordova,  &:  White  (1986)  found  that  the  IR  and  soft 
X-ray  light  curves  had  the  same  phase.  The  consistency 
of  the  hard  X-ray  and  IR  phases  thus  justifies  combining 
the  hard  X-ray  data  with  earlier  soft  X-ray  observations 
to  extend  the  ephemerides. 

Some  caution  needs  to  be  exercised  in  interpreting  this 
result.  Cyg  X-3  is  known  to  have  significant  cycle-to-cycle 
variations  in  the  light  curve  shape.  Here  we  are  compar¬ 
ing  one  IR  cycle  (incompletely  sampled,  and  with  possible 
underlying  variations)  to  the  OSSE  multi-orbit  average. 

5.2.  Light  curve  shape 

Previous  observations  showed  that  the  IR  and  soft  X- 
ray  light  curves  had  the  same  shape  (Mason,  Cordova, 
&:  White  1986).  However,  here  both  the  OSSE  and  IR 
light  curves  are  marginally  inconsistent  (at  about  the  0.2% 
level)  with  the  canonical  soft  X-ray  template  (van  der  Klis 
&:  Bonnet-Bidaud  1989).  Again,  we  cannot  rule  out  that 
the  difference  in  IR  light  curve  shape  was  caused  by  un¬ 
derlying  DC  variations  during  the  measured  cycle  or  by 
the  normal  cycle-to-cycle  differences.  In  fact,  a  short  pe¬ 
riod  of  UKIRT  observations  on  July  16  (one  day  before 
the  light  curve  presented  here)  revealed  Cyg  X-3  to  be  in 
a  significantly  brighter  state  in  the  IR  (Fender  1995).  This 
suggests  a  possible  DC  decline  is  superimposed  upon  the 
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X-ray  template  (dashed  line). 


orbital  light  curve.  This  might  affect  both  the  light  curve 
shape  and  the  fitted  phase. 

Evolution  of  light  curve  shape  with  time  is  expected 
if  the  apparent  orbital  period  change  is  due  to  apsidal 
motion.  Some  evidence  for  this  was  present  in  the  earlier 
IR  observations  of  Jones  et  al.  (1994)  where  the  light  curve 
appeared  more  symmetric,  and  with  a  broader  minimum, 
than  seen  in  previous  measurements. 

The  combined  OSSE  data  offer  some  support  for  the 
hypothesis  that  the  X-ray  light  curve  shape  has  changed 
since  the  1983-1985  EXOSAT  observations.  Figure  3 
shows  the  average  OSSE  light  curve  from  the  seven  ob¬ 
servations  from  1991  to  1994.  This  curve  is  inconsistent 
with  the  template  light  curve,  also  shown  (x2  =  46.8  for 
13  dof,  confidence  level  =  1.0  x  10-5).  The  light  curve  ap¬ 
pears  more  symmetrical,  with  a  somewhat  narrower  min¬ 
imum.  A  transient  change,  over  approximately  one  week, 
of  the  soft  X-ray  light  curve  from  the  template  shape  to  a 
more  symmetrical  form  has  been  observed  by  van  der  Klis 
&  Bonnet-Bidaud  (1981).  However,  such  a  change  should 
have  little  effect  on  the  average  OSSE  light  curve,  which 
combines  ~  65  days  of  observations. 

It  is  possible  that  the  light  curve  shape  has  not  changed 
but  is  intrinsically  different  above  50  keV.  This  cannot  be 
ruled  out  because  there  are  no  comparable  hard  X-ray 
light  curves  from  earlier  epochs.  Also,  there  are  no  cur¬ 
rently  available  soft  X-ray  light  curves  taken  at  the  time 
of  the  OSSE  data  to  make  a  contemporaneous  comparison 
of  light  curve  shapes  in  the  different  energy  ranges.  How¬ 
ever,  the  independent  evidence  for  light  curve  evolution  in 
the  infrared  argues  against  this  alternative. 


5.3.  Modulation 

The  modulation  (light  curve  (max- min) /max)  is  a  key  test 
of  source  models.  In  these  data  we  measured  a  modulation 
of  ~  20%  for  IR  and  ~  50%  for  hard  X-rays.  The  strong 
modulation  in  the  hard  X-rays  is  inconsistent  with  the  ex¬ 
trapolation  of  the  EXOSAT  observations  by  Willingale, 
King,  &  Pounds  (1985),  which  predict  essentially  no  mod¬ 
ulation  at  these  energies.  Strong  hard  X-ray  modulation 
may  be  incompatible  with  their  wind  model  of  the  Cyg 
X-3  light  curve. 

6.  Summary  and  Conclusions 

The  current  observations  show  some  intriguing  possible 
correlations  between  IR  and  hard  X-ray  behavior.  In  par¬ 
ticular,  the  phases  of  flux  minima  are  consistent  at  both 
wavelengths,  and  both  show  possible  deviations  from  the 
historical  soft  X-ray  light  curve  shape.  However,  strong 
variability  in  IR  clouds  the  issue.  To  understand  this 
highly  variable  source  requires  more  coordinated  obser¬ 
vations,  with  better  IR  coverage  over  a  longer  period,  co¬ 
inciding  with  hard  X-ray,  soft  X-ray,  and  radio  measure¬ 
ments. 
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